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A system of a particle and a harmonic oscillator, which have pure point spectra if un- 
coupled, is known to acquire absolutely continuous spectrum when they are coupled by a 
sufficiently strong point interaction. Here the dynamical mechanism underlying this spectral 
phenomenon is exposed. The energy of the oscillator is proven to exponentially diverge in 
time, while the spatial probability distribution of the particle collapses into a <5-function in 
O I- the interaction point. On account of this result, a generalized model with many oscillators 

• which interact with the particle at different points is argued to provide a formal model for 
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approximate measurement of position, and collapse of wave packets. 
PACS numbers: 03.65.Yz: 02.30.Tb 
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A. Background: Smilansky's model. 

Smilansky's model [1] is an offspring of the theory of Quantum Graphs. It consists of a quantum 
particle coupled to a harmonic oscillator via a point interaction. The particle moves inside a 1- 
dimensional hard box . The linear coordinates of the harmonic oscillator and of the particle 
are denoted by q and by x respectively , with x E II = [— L/2,L/2]. In the Hilbert space 
fj = L 2 (Il) ® L 2 (M) the Hamiltonian of the system is formally written as: 

U a ^ = # (p) ®I + \®H^ +aq5(x) , (1) 

where 

1 f) 2 1 

(p) 

is a harmonic oscillator Hamiltonian with frequency u, HY is the Hamiltonian of the particle 
in the box, and the last term describes the point interaction, which is scaled by the parameter 
a > 0. In [1] Smilansky's model was presented in two variants, one with L = +oo and the other 



2 



with L < +00. The spectral theory of the former variant was rigorously analyzed by Solomyak 
[2I and Naboko and Solomyak who proved that for a > oj a new branch of the absolutely 
continuous (ac) spectrum of 7-L a ,oo,u} appears besides the one which is naturally associated with 
unbounded motion of the particle. The new branch coincides with R, and has multiplicity 1. The 
2nd variant is the finite-box model which is studied in the present paper. For this variant an 
argument presented in [I] shows that normalizable eigenfunctions exist for a < oj and don't exist 
for a > oj. 

A generalization of Smilansky's model has N oscillators interacting with the particle at different 
points. Evans and Solomyak [4] have used a scattering theory approach to prove that a spectral 
transition occurs also in the multi-oscillator model with N = 2 and L = 00 . The nature of their 
argument makes it intuitively clear , that the same conclusion is true for any N > 0. 
Not much is known about the dynamics of Smilansky's model. In [ll] strong excitation of the 
oscillator was surmised for a > oj, with the particle dwelling near the interaction point. Due to 
inherent exponential instability of the dynamics, to be proven in the present paper, numerical 
simulation of this system is problematic. 



B. Outline. 



In this paper the dynamics of Smilansky's finite-box model is studied for a > oj, in the single- 
and in the multi-oscillator cases (sees. HT1 IIIII respectively). The main results are Propositions El 
[3] and [5j In the single-oscillator case, the energy of the harmonic oscillator diverges exponentially 
fast, and, in the limit t — > ±00, the probability distribution of the particle collapses into a 
5-function supported in the interaction point. The dynamical origin of such exponential instability 
may be qualitatively illustrated as follows: for q < 0, the q5 term in the Hamiltonian acts like 
a potential well for the particle. When a > oj, interaction drives the oscillator still farther in 
the q < region. This makes the potential well deeper, and so on. Unbounded increase of the 
oscillator's energy follows, which is balanced by unbounded decrease of the energy of the particle 
as it falls deeper and deeper into the well. 

This picture is further illustrated by an approximate description of the dynamics, based on a band 
formalism. This is formally a Born-Oppenheimer approximation in which the particle plays the 
role of the fast degree of freedom; however it is a long-time asymptotic approximation rather than 
an adiabatic one. In this approximation the oscillator gets an effective spring constant, which 
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becomes negative when a > oj. 

The mathematical groundwork for the exact dynamical results of Propositions [2] and [3] is provided 
by spectral results largely resting on the work of Naboko and Solomyak, which are described 
in sections III Al and III Bl In particular, existence of ac spectrum of multiplicity 1 for a > cu 
is assumed as a rigorously proven result, because Naboko and Solomyak's proof ^ of a new 
branch of ac spectrum in the L = oo case works, with minor modifications, also in the finite 
box case. An independent proof of existence of ac spectrum (though not of its simplicity) is 
nevertheless provided here by spectral expansions, which are constructed in section III Bl using 
formal eigenfunctions. Such eigenfunctions are studied in Section III A| by adapting a method 
used in [3], which includes recourse to Birkhoff's theory about asymptotic expansions of solutions 
of 2nd order difference equations In addition, new results about smooth dependence of 

eigenfunctions on energy, which are necessary for the purposes of spectral expansion, are proven 
(Note lVA"]) . elaborating on the formulation by Wong and Li 1J| of Birkhoff's theory. To be noted 
that, unlike the L = oo case, in the finite-box case some point spectrum may survive even above 
the threshold a = u>, in the presence of special symmetries (see section III Aft : and that, even in the 
absence of point spectrum, pure absolute continuity of the spectrum is not proven (nor it was in 
the L = oo case). 

A finite-box model with an arbitrary finite number of oscillators is studied on a somewhat less 
rigorous level. Validity of the scattering approach which was developed in [J] for the L = oo 
variant is assumed, so a spectral transition is again expected and indeed numerical computations 
of bands provide evidence that at a = w the morphology of the lowest bands undergoes a phase 
transition , which mirrors the spectral transition. 

Using the scattering approach, the reduced state of the particle is shown to evolve towards a fully 
incoherent mixture of "position eigenstates" . This process looks like the wave-packet reduction 
which is associated with a measurement of position and indeed the multi-oscillator model is 
surmised to provide a formal model for approximate position measurement, with the oscillators 
acting like detectors of the particle's position 1 . 



1 Models with point interactions , different from those considered in this paper, have already been used in studies 
of decoherence (see, e.g. @, 0|), as well as models of particles interacting with oscillators Q. Approaches to 
decoherence based on scattering have been used e.g. in [y, ||. 
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II. SINGLE OSCILLATOR MODEL. 
A. Formal eigenfunctions. 

For A > let A : ip(x, q) \-t Xip(Xx, Xq) denote the unitary scaling operator from L 2 {Ii) <8> L 2 {W) 
to L 2 (I x -i L ) ® L 2 (R). Then: 

A 1 7i a ,L,u> A = A 2 % a ' t L',u' , (3) 

where a' = a/A 2 , L 1 = LA, and cj' = ujjX 2 . Therefore, one of the parameters L,co,a may always 
be re-set to a prescribed value, by suitably rescaling the coordinates x and q and the time t. Here 
L = 2ir is assumed, and periodic boundary conditions at x = ±ir are used, so the particle may 
be thought to move in a circle § with a distinguished point O. This choice affords some formal 
simplifications without hindering theoretical analysis (see footnote on page 10). That being said, 
a, L, lo will no longer be specified in subscripts to H, unless strictly necessary. As the Hamiltonian 
is invariant under reflection (x \— )■ —x) in the interaction point O , odd functions with respect to 
x make an invariant subspace. Such functions vanish at the interaction point, so in this subspace 
the particle and the oscillator do not interact. For this reason, analysis will be restricted to the 
invariant subspace fj+ = L+(§) ® L 2 (M) where L^_(S) are the square- integrable functions on § 
which are invariant under reflection in O. 

The theory which was developed by Solomyak and Naboko in papers for the case L = oo 

works, with minor modifications, in the present case as well. It starts by analyzing the "formal 
eigenfunctions" of T~L, and this analysis will now be adapted to the present case because these 
very eigenfunctions will provide a key to the dynamical analysis to be presented in the following 
Sections. The wave function ip{x,q) is expanded over the normalized eigenfunctions h n {q) of the 
harmonic oscillator (Hermite functions): 

oo 

ip{x,q) = ^2^n(x) K(q) , (4) 

n=0 

The Hilbert space is thereby identified with ^ 2 (N) (g>L 2 (§), that is the Hilbert space of sequences 
■0 = {ipn(x)} such that \\ip\\ 2 = IIVViH 2 < +°° where |.| denotes the L 2 (S) norm, and the bold- 
face symbol ||.|| denotes the fj-norm. The Hamiltonian is formally identified with the differential 
operator which acts in £ 2 (N) ® L^(S) according to: 

{lpn(x)} l-> {L n tp n (x)} , 

Id 2 ( 1\ 

Ln = ~2d^ + \ n+ 2 r ' ( n = ' 1 ' 2 '--- ) • (5) 
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A vector ip £ ^ 2 (N)®L 2 (§) is in the domain of the Hamiltonian if each ij: n has a square-integrable 
2nd derivative in § \ {O}, such that ^ n ||iy n ^; n || 2 < +00, and, moreover, certain boundary condi- 
tions at x = are satisfied. These are dictated by the 5 function in ([T]). Using recurrence properties 
of the Hermite functions, such "matching conditions" may be written in the form jjj: 

<(0+) -<(0-) = 2a(2u)- 1 / 2 (V^+T^n+i(0) + V^V-n-i(O)) , (6) 

where the rightmost term is for n = 0. The "formal eigenfunctions" are sequences {u n (x,E)} 
that solve the infinite system of equations 

L n u n (x,E) = Eu n (x,E), n = 0, 1, (7) 

for and satisfy the matching conditions ([6]). Such sequences need not belong in £ 2 (N)<g>L 2 (8), 

and are sought in the form: 

u n (x,E) = C(n,E) v n (x,E) , (8) 
where, for each n = 0, 1, 2, . . ., the functions: 

v n (x,E) = p n (E) cos (k n (E)(\x\ -vr)) , k n (E) = yj2E - (2n + l)u (9) 
are normalized solutions of eqn.([7]). The factor p n in © is chosen so that \\v n (., E)\\ = 1 : 

Pn{E) = V + 2k n (E) ) ■ (10) 
For n > E/uj — 1/2, k n (E) is imaginary and formulae ([9]), (|10p are conveniently rewritten with 
circular functions replaced by hyperbolic ones, and k n (E) replaced by Xn(E) = —ik n (E). 
Coefficients C(n,E) in ^ have to be chosen such that conditions © be satisfied. Substituting 
© and © in eqs. © one finds that to this end they must solve the following 2nd order difference 
equation : 

h 2 (n,E) C(n + 2,E) + hx(n,E) C(n + 1,E) + h (n,E) C(n,E) = , (n > 0) , (11) 
with initial conditions C(0, E) and C(1,E) that satisfy: 

M-l,^) C(l,£) = -/n(-l,£)C(0,£), (12) 
having denoted, for n > — 1 : 

^(n,^) = a v n+2 (0,E) ; /n(n,£) = (2w) x / 2 < +1 (0+, £) , (13) 
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and for n > 0: 

h (n,E) = aVn + 1 v n (0,E) . (14) 

Let £ w denote the set of real energies E such that either h2(n,E) or ho(n,E) or both vanish for 
some n > 0: such energies are given by IE = (2n + l)w + (r + 1/2) 2 for r£Z and n > 0, so ^ has 
at most finite intersection with any bounded interval. Whenever E 6 R\£ w , solutions of eqn. (jlip 
are in one-to-one correspondence with their initial values C(0,i?) and C(1,E). In particular, the 
solution of (jlip which verifies (|12p exists, and is uniquely fixed by the value of C(0, E), which plays 
the role of a normalization constant. An asymptotic approximation as n — > +oo to this solution is 



provided by a theory of Birkhoff and Adams 



131 ] , which is briefly reviewed in Note IV Al According 
to that theory, if a > oj and then equation (jlip has two "normal" linearly independent 

solutions C±(n,E) which have the following n — > +oo asymptotics: 

C^(n,E) \= exp(±in# ^i\Elog(n)) + 0(n" 3/2 ) , (15) 



where 

9 = arccos(w/a) , A = — -j^=^= . (16) 

(Cp. Lemma 3.3 in Ref. [3j, where parameters jx and A respectively correspond to u/a and E/uS). 
Thanks to reality of coefficients (113p and (|14p . the complex conjugate of any solution of eqn. (jlip 
is still a solution. Hence the normal solutions are mutually conjugate, because such are their 
asymptotic approximations (]15p . If C(0,E) is chosen real in (]12p . then the sought for particular 
solution C(n, E) is real for all n > 0, so it may be written as a linear superposition of the mutually 
conjugate "normal" solutions, in the following form : 

C(n,E) = C(0,E)$le {Z(E)C\{n,E)} , (17) 

where Z[E) S C has to be chosen so that: 

Me{Z(E)C* + (0,E)} = 1, h 2 (-l,E) Ke {Z(E)CX(1,E)} = -h^-^E) , (18) 



as is required by the initial condition (I12p . In the rest of this work, the normalization factor 
C(0,E) in (|17p will be fixed such that C{0, E)\Z(E)\ = -zr -1 / 2 . This choice is aimed at Lemma 1 
below. Thanks to it, coefficients C(n,E) of the formal eigenfunction {u n (x,E)} (cp. (jSJ) have the 
following asymptotic form: 

C(n,E) ^= cos(n6 - XE\og(n) + ((E)) + 0(n" 3/2 ) . (19) 
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Here C(-^0 is the phase of Z(E); it is not explicitly known, because the normal solutions C±(n,E) 
are not known except for their asymptotic forms, so eqs. (I18p cannot be solved explicitly. In 
Note IVAl (Corollary [2]) Z(E) is proven to be a C 1 function of E in any closed interval I having 
empty intersection with the set £* of "exceptional" energies, which is obtained on adding to £ w 
the threshold energies (n + 1/2)lo, (n > 0), which are branch points for coefficients in eqn. flllh . 

It should be stressed that (fT9j) only holds when a > uj, and that the case a < uj is not discussed 
here. Of crucial importance is that: 

+oo 

Y,\C(n,E)\ 2 = +oo, (20) 

n=0 

so the sequence {u n (x, E)} is not in £ 2 (N) ® L 2 (§>) and does not define an eigenvector proper. It is 
worth noting, however, that the series ^2 n C(n, E)v n {x, E)h n (q) is pointwise convergent to a C°° 
function tp^(x,q) at all points (x,q) G (S \ {0}) x R, because if |x| > 5 > then v n (x,E) decays 
quite fast, ~ 2n 1//4 e _5v/ ^™ with n, and the Hermite functions are uniformly bounded 10]. 
Such properties of the infinite recursion (1111 ) are essentially identical to ones which were established 

n 

in [3( for the L + oo case. Due to them, for a > u the spectrum of % a ,L,u acquires an absolutely 
continuous component, with multiplicity 1. 

B. Spectral expansions. 

Throughout the following, a > uj is understood, and the absolutely continuous subspace of H 
is denoted fj ac . In this Section the above described formal eigenfunctions {u n (x,E)} are used to 
construct spectral expansions. 

Lemma 1 If $ € Cb(K \ (the continuous, compactly supported functions having no excep- 
tional energies in their support), then 

J2 ff dE 1 dE 2 ^!{E 1 )^{E 2 ) [ dxu n (x,E 1 )u n (x,E 2 ) = [ dE V(E)$(E) . (21) 

n-U R2 

Proof, let P n (Ei,E<2) = (u n {E{),u n (E2)) denote the scalar product in L^_(S) of u n {x,E\) and 
u n (x,E 2 ). From L n u n = Eu n (cp. eqn.©): 

P n (Ei,E 2 ) = -(2E 1 )- 1 N + J \dxu n (x,E 2 )-^u n (x,E 1 ) + 

+ (2E 1 )- 1 (2n + l)(u n (E 1 ),u n (E 2 )) . (22) 
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Integration by parts yields: 

2{E 1 - E 2 )P n (E 1 ,E 2 ) = [u' n (0+,E 1 )-u' n {0-,E 1 )]u n (0,E 2 ) + 

- [<(0+, E 2 ) - <(0- E 2 )] u n (0, E x ) . (23) 
Using the matching condition ©, 



{Ex - E 2 )P n {E 1 ,E 2 ) = — (W n+ i(Ei,£ 2 ) - W n (^i,S2)) , (24) 



where, for n > 

W n (^i,^) = ^K(0,£iK-i(0,£ 2 ) -Un-iCO^iKCO,^)) (25) 
and W (-Ei,-Ba) = 0. Hence, 

AT 

rP n(£l ,E 2 ) = *W N+1 (E U E 2 ) 

Substituting in (j25[) the asymptotic form of u n (Q) which follows from eqs,(f9]). ()10p . (|19p : 

^ p ,„ 9 . ... sin((A(E 1 - £ 2 ) ln(AQ - CgQ + C(^)) 
2^n(£l,£ 2 ) ~ 2a sm(6l) WETT^ 

n=0 

S(E X -E 2 ). (26) 

N— >oo 

where the definitions (|16p of A and 6* have been used. Convergence to the Dirac delta function is 
meant in the sense of eqn. (|2ip for continuous functions and $ supported in R\£* and rests on 
regularity of ((E) in R\ £ *, as established by Corollary [2] in Note IVAl D 

Thanks to (|2ip . whenever ^ E Cq(R\£*) the sequence of functions which are defined on § by 

Mx) = [ dE*(E)u n (x,E) , (n = 0,1,2,...) (27) 

is a vector in £ 2 (N) ® Lj_(S). This vector will be denoted ip, and the function will be termed the 
spectral representative of ip. Eqn. (|2ip says that the map i : i— » ^ is isometric, so z extends to 
an isometry of L 2 (M) into £ 2 (N) ® £+(S) . Proposition Q] below easily follows . It shows that this 
map, or rather its inverse, yields a complete spectral representation of H restricted to its absolutely 
continuous subspace. 

Proposition 1 ; 

(i) For all V E L 2 (R), and t E R, 



~ iEt q>) = e - iW i(ty) , (28) 



(ii) i is a unitary isomorphism of L (R) onto Sj ac . 



A proof is presented in Note IV Bl 
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C. Dynamics for a > uj. 

The spectral results in the previous Sections have dynamical consequences as stated in the 
Propositions below. Throughout this section a > uj is understood. Let tp G <8> L 2 (M.) , and 

11-011 = 1. The notation ip(t) = e'^ip will be used; moreover, {/)t = y fn ' dtf(t) will denote the 
time average up to time T of a function f(t). 

Proposition 2 If ip € Sj ac then the time-averaged energy of the oscillator grows in time, at least 
exponentially fast: i.e., 

\n(E osc (T)) 

hmmf — > 0. (29) 

where 

E osc (T) = ( {m , I®H^i>(t)) 

\ I T 

A proof is given in Appendix IV CI 

Proposition 3 If ij) £ Sj ac and ||^|| = 1 then the probability distribution of the position x of the 
particle weakly converges to 5(x) in the limit t — > +oo. 

This is equivalent to 

lim dq / dx \ip(q, x, t)\ 2 = (30) 

t^+co J J 

r/<\x\<n 

for any < rj < tt, which is proven in Appendix IV Dl 

One may reasonably expect the expectation value of the position q of the oscillator to diverge to 
—oo as the particle endlessly falls in the ^-potential . This will be further supported by arguments 
in the next Section; however no exact proof is attempted here. 

The reduced state of the particle when the full system is in the pure state tjj(t) is the positive trace 
class operator S^(t) in L 2 (S) such that Tr(S^(t)A) = (rJ)(t),A(gihJ)(t)) for all bounded operators 
A in L 2 (§>). Proposition [3] entails a somewhat extreme form of decoherence for the reduced state: 

Corollary 1 If tp G f) ac then for every <p and <p' in L 2 (S), lim t ^ > .± o((f), S^,{t)(j)') = . 

A proof is presented in Note IV El 
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FIG. 1: The lowest energy band in the single oscillator model for u> = 1 and a = 0.8 (upper, full line), a = 1 
(middle, dotted line) and a = 1.3 (lower, dashed-dotted line). 

D. Band dynamics: inverted oscillator. 



An intuitive picture of the above results is provided by an approximate description of the 
dynamics, to be presented in this section. Hamiltonian (pQ) (with L = 2tt) may also be presented 
in the following form: 



where, for any fixed value of q, 



H = / dqB a (q) + I (g) H^ sc> 
JR 



H a ( 9 ) = + aq5{2 



(31) 



(32) 



is an operator in [ill ]. It has a complete set of eigenfunctions and eigenvalues which para- 

metrically depend on the product aq. All eigenfunctions are real valued and have the form: 

(f>g,n(x) = A n COs(i n (\x\ - 7r)) 71 = 0,1,2,... (33) 

where A n are normalization constants, and £ n are the solutions with 9?e(£) > of the equation: 

tan«) = — . 

? 

They are numbered in increasing order of the corresponding energy eigenvalues W n (q) = ^Cn- ^ n 
£ n with n > are real, and behave like £ n ~ n + | asymptotically as q — > ±oo. Instead is rea l 
only when q > 0, and turns imaginary when q < ; in that case £o = i\, where \ ls the unique 
positive solution of 

tanh(7rx) = — — . 

X 
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Hence Wq (q) is negative whenever q < 0. A standard Feynman-Hellman argument yields: 

dW n (q) 2 

— ^ — = a<£ g ,n(0) , (34) 

so the levels W n (q) are nondecreasing functions of q. The ground state energy Wo(q) is asymptot- 
ically given by: 

s 

1 



W (q) = i£o ~ < 



for q — > +oo ; 

(35) 



\a 2 q 2 , for g — > -co . 



2' 

When q > 0, the ground state eigenfunction has still the form (|33p with n = 0. For q < it is 
instead given by: 

4> g ,o(x) = A cosh(x(<?)(|x| -7t)) , (q < 0) . 
so for large negative g it is sharply peaked at x = 0. Any ^ 6 £+(§) 8) -L 2 (M) may be expanded as 
■0(ar,g) = ^ Qn(<?) ^ 9 ,n(^) , Qn{q) = / dx (j) q>n (x)i/j(x,q) , (36) 

n=0 •^ 7r 

so that 

OO /*7T 

|2 



ElQn^l 2 = [* dxMx,q)\' 

n=0 •^ 7r 



In this way fj is decomposed in "Band Subspaces" 5B n = {VK^j ?) = Qn(q)4>q,n(x) \ Q n (q) £ L (M)}. 
The projector onto the n-th band subspace will be denoted II n . In the "band formalism" it is easy 
to show that T-L a ,u is bounded from below when a < u. Indeed, if ip is in the domain of 7i a ,w, then 

°° /"OO 

{^H a ^) = / dqW n {q)\Q n {q)\ 2 + (^>, I ® fl^V) • (37) 
Singling out the contribution of the lowest band 23o; and using ([35]) and monotonicity of Wo(q): 
{i/>,'Ha,ui>) > -^« 2 / dq I dx q 2 \tp{x, q)\ 2 + 



°° /"OO 

+ J2 d Q Wn(q)\Qn(q)\ 2 + tf<T c) lf>) 

n=l- / -°° 



> ly/ u 2- a 2 ||^f . (38) 

This shows that the abrupt change from semibounded to unbounded spectrum which occurs at 
a = uj is related to a change in the structure of the "ground band" alone. This transition is further 
elucidated by noting that the ground band plays the role of a stable variety, due to the following: 
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Proposition 4 If a > ui then, for arbitrary ip S Sj ac; and for all integer n > 0: 

lim || Jl n e~ iHa ' ut ip II = . 



t— >oo 



A proof is presented in section IV Fl This suggests that asymptotic solutions in time of the time- 
dependent Schrodinger equation may be sought in the form ijj(x,q,t) = Qo(q,t)4> q fi(x). No exact 
proof is attempted here; nevertheless, direct substitution yields a Schrodinger equation for the 
band wavefunction Qo(q,t): 

MQo ld 2 Q 



dt 2 dq 2 

where the band potential V(q) is given by: 



+ V(q)Q Q 



VUl) - \" 2 q 2 + W (q) + \J dx 



(39) 



dq 

A calculation reported in Note IV Gl shows that the last term on the right-hand side is 0{q~ 2 ) as 
q — > — oo and tends to a constant when q — > +oo. At large q, the band potential is then determined 
by the other two terms. The band potential which results of these two terms alone is shown in 
Fig.l. As a grows beyond uj, it turns from concave to monotone increasing and then, at large 
negative q, V(q) ~ — \[c? — uj 2 )q 2 , that is the potential of an inverted harmonic oscillator. This 
sort of phase transition qualitatively explains the growth of the harmonic oscillator's energy which 
was proven in Proposition 1, and suggests that it is exponential with rate A" 1 = 2\/a 2 — lo 2 . 

III. COLLAPSE OF WAVE-PACKETS. 
A. A multi-oscillator model. 

Smilansky's model has generalizations, in which an arbitrary finite number of harmonic 
oscillators interact with the particle at different points Oi, . . . ,On ■ In the L = oo case, the 
corresponding spectral theory has been developed by Evans and Solomyak 4|] using a scattering 
theory approach. Translated to the present case, this approach is as follows. All oscillators are 
assumed to have the same frequency u and coupling constant a, and a circular ordering is assumed 
for the interaction points Oi. For each i = 1, . . . , TV let a rigid wall be inserted at a point Z{ in 
between Oi ad Oj+i, and let Ji denote the arc The Hamiltonian %W of the resulting 

system differs from % because of Dirichlet conditions at the points Zi, and is actually an orthogonal 
sum of operators np in S)i = L 2 (Ji) 8) L 2 (M. N ), each of which describes the particle in a rigid 
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box J;, coupled to the i-th oscillator alone. Therefore, thanks to what is known about the single- 
oscillator box model, T-L^ has ac spectrum coinciding with R when a > uj. M0ller wave operators 
are defined by : 

n±(H,H^) = lim e^e-^'pW (40) 

t— >±oo 

where denotes projection onto the absolutely continuous subspace of T-L^ h \ They are said to be 
complete if their range coincides with the entire absolutely continuous subspace of T~L. Whenever 
this happens, the wave operators U±(T-C^ b \T-C) also exist [ijj]. Existence and completeness have 



been proven in 



4, 



for the model with L = oo. Here they are assumed also for the multi-oscillator 
box model. In the absence of a formal proof paraphrasing Evans and Solomyak's, this assump- 
tion rests on intuition provided by Proposition 2: as the wave function which evolves in fjj under 
Hamiltonian is drained by the i-th interaction point, the boundaries at Z{ and Zi + \ become 
ininfluent, and so does the difference between T~L and ? 

Existence and completeness of wave operators enforce unitary equivalence of the absolutely con- 
tinuous parts of H and of 7i^ b \ hence infinitely degenerate Lebesgue spectrum of % at a > uj. 

B. Band formalism. 

There is a band formalism also for the N-oscillators case. The case N = 2 will be briefly 
described. Oscillators 1 and 2 with respective coordinates q\ and qi are coupled to the particle 
at points x = and x = ±7r respectively, diametrally opposite in §. The particle Hamiltonian 
which now replaces (|32]) parametrically depends on q\ and q2, and has real- valued eigenfunctions 



<f>q lt q 2 ,n( x ) m -^ 2 (§>+) an d eigenvalues W n ((71,(72) = f£n> where £ n are the solutions with 3fte(£ n ) > 
(numbered in non-decreasing order of the corresponding eigenvalues) ) of the equation: 

tan« = — . (41) 

i l - a z qiq 2 

"Band potentials" E n (q\,q2) = \u 2 {q\ + q%) + W n (qx,q2) computed by numerically solving eq. 



(I4ip are shown in Fig. 2 and in Fig. 3. Like in the N = 1 case, the spectral transition at a = uj 
is concomitant to a phase transition in the structure of the lowest band. A transition is observed 
for the 2nd lowest energy band as well, at a higher value of a/ui ~ 1.414 (not shown), but not 
for higher bands. The structure of the overcritical ground band, shown in Fig. 3, is explained as 



2 The same picture accounts for irrelevance of boundary conditions in the single-oscillator box model. 
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follows. When 1 < a/uo , £o is found to be imaginary in the region 1Z of the (91, q 2 ) plane which is 
defined by the inequality: 



where q- and 9+ are respectively the minimum and the maximum of q\ and q 2 . Moving out to 
infinity in 1Z along a half-line started at (0,0), the asymptotic behaviour of £0 is ~ —ictq—, so 
the band potential Eo(qi,q 2 ) diverges to —00, except along directions lying within an angle of 
arcsin(o;/a) — it/ 4 on either side of the half-line qi = q2 < 0, where instead it diverges to +00 as 
long asw<a< ojy/2. This is why in Fig. 3 one observes two valleys , hereby labeled 1 and 2, that 
descend to —00 along the negative q\ axis and the negative 92 axis respectively, and are separated 
by a crest, which rises along the q\ = q 2 < half-line. In region 1Z the ground eigenfunction has 
the form 

Kmfl{x) = C x {q x ,q 2 ) e^>v)\W + C7 2 (?i, q 2 ) e -^^ x \ . 

It has two peaks, labeled 1 and 2, at the interaction points of oscillators 1 and 2. Descending along 
either valley both peaks become narrower and narrower, however calculation shows that the whole 
probability is asymptotically in time caught within the peak which shares the label of the valley. 
Like in sec. Ill Dl for oj < a < wV2, the quantum dynamics is asymptotically attracted by the 
0-th band subspace, which consists of functions of the form Qo(qi, q 2 )<fio,q lt q 2 (x) ■ The "band wave- 
function" Qo(q±, q 2 , t) asymptotically in time solves the Schrodinger equation for a particle in the 
plane (gi, q 2 ), subject to a potential that behaves at 00 like the one in Fig. 3. A classical particle 
would escape to infinity along one valley, so one and just one oscillator would undergo unbounded 
excitation. 



C. Collapse of wave-packets. 

Existence and completeness of wave operators have the following immediate consequence, which 
generalizes Proposition [3l 

Proposition 5 If a > uj and t/j G % a c with ||^|| = 1, then the probability distribution of the 
particle converges weakly as t — >■ ±00 to a superposition of 5 functions supported in the interaction 
points. 

j j ' dqi...dq N \il){x,qi,...,q N ,t)\ 2 ^ -yf 5(x - Oj) , 
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FIG. 2: The lowest energy band in the model with 2 oscillators, for a = 0.7 and u = 1. 




-100 _ioo 



where: 



FIG. 3: Same as FigH for a = 1.3 and w = 1. 



7, 



|/>-ft ± (^ 6 \W|| 



(42) 



and Pj denotes projection onto Sjj . 

A proof is presented in section IV HI This in particular implies that the right-hand sides in ()42|) do 
not depend on the positions Zi of the rigid walls. 

Corollary [1] about complete decoherence of the reduced state of the particle generalizes to the multi- 
oscillator case. Hence one may say that coupling to the oscillators causes the reduced state of the 
particle to evolve exponentially fast towards an "incoherent superposition of position eigenstates" . 
In the case when the particle is initially in a pure state, this process is similar to the wave- 
packet reduction which is conventionally associated with measurements of position. Here the 
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measuring apparatus consists of N oscillators, and the iV interaction points have to be chosen in 
a thick homogeneous grid. Under the assumptions in Proposition [5j as t — > +00 the pure state 
comes closer and closer to the state e l ^ mt Q + ifj. This is a coherent superposition 
of states which have the particle in a box and the corresponding oscillator in a highly excited state. 
Tracing out the oscillators yields an incoherent mixture of alternatives for the particle position , 
however the probability of finding the particle in the m-th box is given by 7+ = ||-P m f2+ , 0|| 2 as 
in eqn. (f22j) and not by fj dx\ipo(x)\ 2 = \\P m iJj\\ 2 as in an ideal measurement. The difference 
lies with replacing the initial state tp with the "outgoing state" £l + ip , and may be ascribed to 
the non-instantaneous nature of the measurement process. Increasing the number N of oscillators 
(hence increasing the precision of the measurement) while keeping a/N and uj/N constant causes 
the time scale of the exponentially fast reduction process to decrease proportional to 1/N (see the 
scaling rule and remarks in the end of sect ill Dl) . This suggests a possibility of retrieving the 
ideal measurement of position by a suitable limit process. 

IV. CONCLUDING REMARKS. 

Dynamical instability in Smilansky's model is due to a positive feedback loop between fall of 
the particle in the 6- potential well and excitation of the oscillator. This effect may not crucially 
rest on point interaction, nor on linear dependence of the interaction on the coordinate of the 
oscillator. While such special features are probably optimal in simplifying mathematical analysis, 
a similar behavior may be reproducible with smoother interaction potentials and also in purely 
classical models. 

Smilansky's model is somewhat unrealistic from a physical viewpoint, as it is not easy to conceive 
of physical realizations, albeit approximate. Generalizations of the model to higher dimension, and 
more realistic couplings - if at all possible - may enhance physical interest. 

Acknowledgment: I thank Uzy Smilansky for discussions about his model and Raffaele Carlone 
for making me aware of exact results in related fields. 
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V. NOTES, AND PROOFS. 
A. BAfcWL Theory. 

Corollary [2] to Proposition [U which is proven in this Note, is an essential ingredient in the 
derivation of the spectral expansion in Sect. Ill Bl (notably in the proof of Lemma [1]). Proposition 



1J| in the 



[6] is proven by a rephrased version of a method which was introduced by Wong and Li 
context of a theory of Birkhoff and Adams about asymptotic expansions for n — > +oo of solutions 
of 2nd order difference equations of the form : 

C(n + 2) + p(n) C(n + 1) + q(n) C(n) = . (43) 

The main result of that theory (Theorem 8.36 in ref.[13(]) is that , whenever coefficients p(n) and 
q{n) have asymptotic expansions for n — > oo in powers of n : 

+oo +oo 

p( n ) ~ ^2 a ( k ) n ~ k > q( n ) ~ 5Z&(*0 n ~ k > ( 44 ) 

fc=0 fc=0 

the equation has two linearly independent "normal" solutions, which have asymptotic expansions: 

oo 

C^in) ~ a^n Q± ^ c± (s)n- s (45) 

where a± are the (assumedly distinct) roots of the equation: 

a 2 + o(0) a + 6(0) = , (46) 

and 

0(1) q ± + 6(1) 
± a (0) a± + 26(0) ' 1 ' 

Coefficients c±(s) in (f4"5j) are recursively determined by directly substituting ([43|) in (j33]) with 
c ±(0) = 1- Eqn. (|lip may be written in the form of eqn. (|43|) . with coefficients that additionally 
depend on E: 

p(n) = pin,E) = - h -^± , g (n) = q (n,E) = - , (48) 

h 2 in,E) h 2 in,E) 



where ho, hi, h 2 are as in eqs. (|13|) and (|T4|) . Using ([9]) and (|T0|) one computes asymptotic 
expansions (|44p . In particular, 

a(O,£0 = 17, o(l,^0 = + KO,£0 = !> KhE) = -1, 
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whence it follows that expansions (|45p have the form (|15p at lowest orders. 

In the following / will denote an arbitrary closed interval contained in M \ £ *; positive quantities 
only dependent on a, u), I will be denoted by cj, C2, . . .; derivatives with respect to 2£ will be denoted 
by a dot, like, e.g., in p(n,E),q(n,E) . . .. The following Lemma [2] sets premises for the proof of 
Proposition [6j 

Lemma 2 For all positive n, p(n, E) and q(n, E) as given by M'J\) . j 1$ are C 1 functions 

of E 6 /. Their asymptotic expansions |.^[ ) are uniform in I . Their derivatives have uniform 
asymptotic expansions in I in powers of n _1 , with coefficients given by the derivatives of the 
coefficients a(s,E), b(s,E), as specified in f^<g| j ; &13\) . and (T^Jj. In particular, n\p(n,E)\ and 
n\q(n,E)\ are bounded in I by some c\ > 0. 

Proof, by direct inspection. □ 

Proposition 6 For all n > the normal solutions C±{n,E) of eqn. ft43\), with coefficients as in 
{4^ > are C l functions of E £ I. 

Proof, the proof is the same for both normal solutions, so suffixes ± will be left understood 
throughout. Thanks to Lemma [2j all coefficients c(s,E) are C 1 functions of E £ /, because each 
of them is determined by a finite number of coefficients a and b. Let 9\ and 9\q denote operators 
that act on sequences w : N x R — > C according to 

{mw)(n,E) = w{n + 2,E) + p(n, E)w(n + 1, E) + q(n,E)w{n,E) , (49) 
(Mow)(n,E) = w(n + 2,E) + a(0, E)w(n + 1, E) + 6(0, E)w(n, E) . (50) 

Let Y(n, E) = exp(±m(9 =F i\Elog(n)), and let a normal solution be written in the form: 

C*(n,E)=L N (n,E) + e N (n,E) , (51) 

where ./V is an integer, and L^(n, E) is obtained on truncating at the (N — l)-th order the asymp- 
totic expansion (I45P of the normal solution: 

JV-l 

L N (n,E) = n- 1/2 Y(n,E) c(s,E) n - s . 

o 

Direct calculation yields 

(fRL N ){n,E) = n- 1/2 Y(n,E)R N {n,E), (52) 
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where , for any fixed n, Rjsr(n, E) is a C 1 function of E G /, because such are all coefficients c(s, E), 
(0 < s < N — 1); and, moreover, 

R N (n,E) = 0(n~ N - 1 ) , (53) 

uniformly with respect to E £ I as n — > oo. Differentiating (|52p on both sides, Rn(ti, E) is found 
to have a uniform asymptotic expansion in / in powers of n _1 , so (|53p entails that 

iW(n,£) =0(n- JV - 1 ) , (54) 

uniformly in /. Substitution of ([5T]) and ([52]) into (jl3|) yields: 

(^ eiV )(n,S) = -n" 1/2 F(n,^) R N (n,E) , 

which is equivalent to 

(<Roejv)(n,£) = -n~ 1/2 Y(n, E) R N (n, E) - q{n, E)e N (n, E) - p(n,E)e N (n + l,E) , (55) 

where p(n,E) = p(n,E) — ao(E) and q(n,E) = q(n,E) — bo(E) . Eqn. (j55|) may be read as a 
inhomogeneous 2nd order difference equation, so it can be rewritten in "integral" form using a 
"Green function" for the operator D^o- This is provided by the function 

, sin((n- 1)9) , . 

g W =.(»-!) " Mr) >> (56) 

where s(n) = 1 for n > 0, and s(n) = for n < 0. Therefore, introducing the operator © that 
formally acts on sequences as in: 

+oo 

(<8w)(n, E) = - ^G{n- k) [q(k, E) w(k, E) + p{k, E) w{k + 1, E)) , (57) 

k=n 

and denoting 

+oo 

d N (n,E) = -^G{n-k) k^/ 2 Y(k, E)R N (k, E) , (58) 

k=n 

the sequence e]\r(n,E) must solve the following equation written in vector form : 

e N = d N + e N . (59) 

No solution of the homogeneous equation 9^o£Ar = appears on the rhs of (f59|) . because such 
solutions do not vanish at infinity, as is instead required of e]y(n,E). Let a and denote the left 
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shift operator and its adjoint: (aw)(n,E) = w(n + 1,E), (n > 1), (a'w)(n, E) = w(n — 1,E) if 
n > 1, and (a^w)^, E) = 0. If ejq satisfies (|59l) . then £n := o~ N satisfies: 

e N = a N d N + e N e N , <5 N = a N ®^ N . (60) 

The operator &n is explicitly given by eqn. (i57]) after replacing a, b by a N a, a N b respectively. 
Thanks to Lemma [3] and to the Contraction Mapping theorem, if N is sufficiently large then 
eqn. (|60l) has a unique solution in the Banach space Xn,i of sequences w : N — >• C l {I) such that 

'■= \\ w \\n + ||w||jv* < +oo , (61) 

where 

|H|jv = sup {(N + n) N+1/2 \w(n, E)\ , n > 0, E e 1} , 

\\w\\ N . = sup l(N + n) N+1 / 2 - | M n,E)\ , n > 0, Eel) . (62) 
[ log(iV + n) J 

The thus found £n determines €n{ti,E) as a C 1 function, and hence, via eqn. (!5ip . the normal 
solution, for n > N + 1. For such n the thesis is then proven, because L]\r(n,E) is itself C 1 wrt E 
thanks to already noted properties of coefficients c(s,E). The values of the normal solution thus 
found at n = N + 1 and n = N + 2 can then be used to retrieve the normal solution for < n < N 
by solving eqn. (143p backwards (which is possible, because q(n,E) ^ for all n > 0, thanks to the 
assumption that E is not in £ *). As this process involves a finite number of steps, and p, q are C 1 
functions, the proof is complete. □ 



Lemma 3 (i) o djy G %n,i, (ii) ©at is a bounded operator in Xn,i, an d its norm is bounded by: 

||0|| £ < 2ec 2 (ci + 3/3)(iV + l/2)- 1 , 
where c\ is as in Lemma\^ C2 = sin(#) -1 = (1 — uj 2 / a 2 ) -1 / 2 , anrf 

/3 = sup {k (\p(k, E)\ + \q(k, E)\ , fe £N, E G /} . 
Proof: (i) from eqs. ([58l) and (j55]l : 

+oo 

div(n,^) < ^^A;-^ 3 / 2 = 0(n- N -^ 2 ), 

k=n 

therefore 1 1 c^iv 1 1 jv is finite and so is ||c"' v+1 djv||Ar. Next, the derivative of the k-th term in the sum 
on the rhs in eqn. ([58]) is 0(A;- JV - 3 / 2 log(l + k)), so the sum of such derivatives is absolutely and 
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uniformly convergent in I to the derivative of d]\r(n,E), and ||<i/v||jv* is finite, 
(ii): noting that 

sup{(£; + N)a N p(k, E) \ k > n} < p , 
and similarly for q, one may write: 

+00 

\(<8 N w))(n,E)\ < c 2 (3\\w\\ N Y,(N + k)- N - 3 / 2 

k=n 

r+00 

< c 2 /3\\w\\ N / (N + xr N - 3 / 2 
= c 2 f3\\w\\ N (N + l^)- 1 ^ + n - ^ (63 ) 

so 

\\& N w\\ N < 2c 2 e/3(iV + l/2)- 1 ||u'|| i v (64) 

thanks to (N + n) N+l / 2 {N + n- i)-^-V2 < 2e. To estimate \\(<8 N w)\\ m : 

+00 

(& N w)\n,E) < caCil^llTv^^ + iV)"^ 3 / 2 + 

k=n 

+00 

+ c 2 p\\w\\ N * Y,( k + N)- N -V 2 log(k + N) . (65) 

k=n 

Estimating the sums on the rhs as it was done in (I63p leads to: 

|| (<3 N w)\\ jv* < (2dc 2 e||ti;||j V + 4c 2 e/3|H|jv*) (iV + 1/2)" 1 . 

Thanks to definition ()61 j) , the latter estimate along with (|64|) yield the claimed bound on the norm 
of ©at as an operator in Xn,i- D 

Corollary 2 If E £ i/ien i/ie difference equation has a particular solution which satisfies 

the initial condition ( tii?j) . and moreover has the asymptotics I119\) . where C(E) *s a C 1 function of 
E in any closed interval of energies containing no exceptional points. 

Proof: the complex amplitude Z(E) (cp. eqn. (|17p ). which determines the sought for solution in 
terms of the normal solutions, is found by solving eqs. ()18p ; so it is a smooth function of the values 
of the normal solutions at n = and n = 1. The conclusion follows, because C(E) is the phase of 
Z{E). □ 
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B. Proof of Lemma [TJ 

First it will be proven that if E m ^(E) G L 2 (R) for some integer m then ip = is in the 
domain of H m , and i(E m ^) = H m ip. 
It is easy to see that 

I%i/j n (x)=c [ dE-$(E)L™u n (x,E) = c [ dE E m ^(E)u n (x,E) 
Jr Jr 

(with L n defined as in ([5])) holds for all Vl/ G Cq(M \ £*) and all positive integers m,n; so, 

thanks to ([21]) the sequence {L™xp n } n is in £ 2 (N) ® L%(&) whenever E m ^(E) G L 2 (R). On the 

other hand the sequence {L™u n } n satisfies the matching condition ([6]) because so does u n , and 

because L n u n = Eu n . The same is then true of the sequence {L™^ n } n , because the condition 

Vl/ G Cb(M\£*) allows for computing left- and right-hand derivatives of L rn "ip n (x) at x = under the 

integral sign . Therefore ip is in the domain of T~L m whenever VP G Cq(R\£*), and 7-L m, if) = i(E m ^). 

As H m is a closed operator, the same is true whenever \& G L 2 (R) and E m ^ G L 2 (R). 

(ii) follows by continuity, because % is isometric. 

(hi) To prove that i is onto: thanks to (|21 j) and (|28p . the time-correlation (ip,e~ l ^ t ip) coincides 
with the Fourier transform of ^(-E)! 2 . Therefore, ^(E)! 2 is the density of the absolutely 
continuous spectral measure of ip with respect to % (also known as the local density of states). As 
% has a simple absolutely continuous spectrum coinciding with R, ip is a cyclic vector whenever 
its local density of states is Lebesgue- almost everywhere different from zero. So, whenever ^(E) 
is a.e. nonzero, ip is a cyclic vector of %, so the closed span of {e~ Mt ip} t£ ^ is the whole of 
£ 2 (N) ® L^.(S), whence ^ 2 (N) ® = *(£ 2 (R)) follows. 



C. Proof of Proposition [2j 

The expectation value of the energy of the oscillator in a state ip = {ip n (x)} G £ 2 (N) <8> •£+(§) of 
the composite system is given by: 



n=0 

and so the sequence 



OO /»7T 

foM®#<r^) = £>+|)w / dx|^(x)| 2 , 



/7T 
dx |u n (x,£)| 2 , (n>0), (66) 
-7T 
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may be thought of as a non-normalizable distribution of the energy of the oscillator over its un- 
perturbed levels, when the full system has energy E. The present proof of Proposition [2] rests on 
the following inequality, which is an immediate consequence of eqs. (j8l).(j9l), (fTOl) . and (fl~9l) : 



limsup nP(n,E) < tt 1 . 



(67) 



Let ip G f) ac , ||V>|| = 1, and ^(E) its spectral representative, so that ^(S)! 2 is the density of 
the absolutely continuous spectral measure of ip with respect to H. Thanks to ([67"]) there is a 
continuous function ^i(E) , compactly supported in M\ £*, so that on the one hand: 



J dE \V(E) -^i(£)| 2 < 1/8 , 



and on the other hand : 



N 



P(n, E) <C ln(iV) , ViV £ N 



(68) 



(69) 



n=0 



for some positive constant C and for all E in the support of ^>\. Let ip\ = x(^E f i ) , ip2 = ip — ipi, and 
^(t) = e~ mt i) = {ip n (x,t)} £ £ 2 (N) ® L 2 + (S). For T > , the probability of finding the energy of 
the oscillator in its ra-th level, averaged from time to time T, is: 



Pn(T) = - / ill 



dx \tp n {x,t)Y 



(70) 



Let pi >n (T) and p2,n(T) denote the functions which are defined by the same equation, with ip 
replaced by ipi and ip2 respectively. By construction of Vl/i, the spectral representation of ipi has 
the form ([2T]) . so Proposition [1] yields: 



i J dt j dx J dEe- iEt ^i(E)u n (x,E) 



(71) 



Vl,n{T) = 7f 

On the other hand, denoting J 7 ^] the Fourier-Plancherel transform in l?{ 

/ dt dE e~ iEt y>x{E) u n (x,E) = 2ir dt u n (x, .)] (t)\ 2 (72) 

Jo J Jo 

< 2vr||^[^ lUn (x,.)]|| 2 = 2ir||*i u n (x, .)|| 2 (73) 
= 2tt J dE ul(x,E)\*i(E)\ 2 . (74) 

Replacing (|74l) in ([71]) . and using ([66]) and inequality ([69]) . which holds throughout the support of 

N 



5>l,nC0 < 2vrCT- 1 ln(iV) . 



(75) 



n=0 
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With N = N(T) := e ClT , where d = l/(16vrC), this estimate yields: 

E P^(T) < h (76) 

n<N T 

From the definitions of p n (T),pi^ n (T),p2, n (T), and ineq. (18) it immediately follows that: 

JV(T) N(T) N(T) 

E ^( T ) ^ 2 E P^ T ) + 2 E ^,n(T) < i + i = i , 

n=0 n=0 n=0 

and so 

i f dt^^H^m) > (l + e > W :z > = i eCiT ^- ( 77 ) 

,/0 n>N(T) 

D. Proof of Proposition [3j 

For < 77 < 7r let = (— vr, — 77) U (77, 7r). The explicit form of u n given in eqs.([8]) and Q 
shows that |u n (a;)| < \C(n, E)\p n cosh.(x n (E)(-K — 77)), whenever n > E/co — 1/2 and x E A7; so, if 
in addition then from (jXOj) and from the asymptotic formula (|19p it follows that 

+00 „ 

E / dxu 2 n (x,E) < +00. (78) 

J A v 

because the integrals in the sum decrease exponentially fast for n — > 00. Thanks to (|78p . for 
< e < 1 one can find a compact set B e C R \ £ * and a continuous function ty e supported in B t , 
so that, on the one hand: 

+00 „ 

E/ dxu 2 n (x,E) <C e , ME E B e , (79) 

for some positive constant C e ; and, on the other hand, — ifi^W 2 < e, where ip^ = i(^f e )- Then: 

dxdq mx,q,t)\ 2 < 2/7 \i>(x,q,t)-^(x,q,t)\ 2 + 



+ 2 jj dxdq \ifj {e) (x,q,t)\ 2 
i,xM 

< 2||^ -V (e) H 2 + 2 ffdxdq \^ e \x,q,t)f 



< 2e + 2 // cfcdg h/) (e) (a;,g,t)| 2 . (80) 
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From eqn.(U]): 



// dxdq\^\x,q,t)\ 2 =J2 dx\iPt\x,t)\ 



and then, since the spectral representative \E' e of ip( € > is compactly supported in R\£*, the spectral 
representation (j27|) can be used to the effect that: 



JJdxdq \^(x,q,t)\ 2 = jj dE dE' e^'"^* V e (E)V e (E') G V (E,E') , (81) 



B £ xB E 

where : 

+00 



u n (x, E)u n (x, E') . 



lr l n=0 

Thanks to (|79|) . G^E^E') is bounded in B e x P e ; on the other hand *$> e (E) is summable over B t , 
so the integral on the rhs in (|8ip tends to in the limit t — > 00 thanks to the Riemann-Lebesgue 
lemma. From (1801) it follows that 



limsup // dx dq \\fj(x,q,t)\ 2 < 2e , 
t->oo J J 
A„xR 



whence the claim (|30p follows, because e > is arbitrary. 

E. Proof of Corollary QJ 

By positivity of S(t), it is sufficient to prove the claim for = 0' and ||0|| = 1. Let P^ = 
(<ft, .)4> denote projection along 4>, and let P v , P^ respectively denote projection onto the functions 
supported in r] < \x\ < tt, and its orthogonal complement. Then: 

(<l>,S(t)<j>) = Tr(5(t)P ) < |Tr(5(t)P P,)| + \Tr (S(t)P+Pf)\ . (82) 

For any r\ > the 1st term on the rhs in ([82]) tends to as t — > ±00 thanks to eqn. ([30j) . due to 

|Tr(S(i)P P,)| = \(ip(t),Pj,P ri ®I'ip(t))\ < ||P„®I^(t)|| . 

On the other hand, the 2nd term on the rhs in (j82|) can be made arbitrarily small, uniformly with 
respect to t, by choosing 7/ small enough: 

|Tr(S(i)P^)| < ||P^|| < ||P^|| , 

Hence the lhs in (f82|) . which does not depend on rj, tends to in the limit t — > ±00. 



26 



F. Proof of Proposition^ 



\Qn(q,t)\' 



dq 



\x\<rj Jri<\x\<TT, 



< 



dq{R( V ,q) + S( V ,q)} 



\x\<7) 



dx (/> q , n (x)ip(x,q,t) 



ri<\x\<w 



dx (f>q 7 n(x)il>(x,q,t) 



where: 

R(r),q) = 2 

s(v,q) = 2 

From the Cauchy-Schwarz inequality: 
/ dqR(ri,q) < 2 [ dq 

JR Jr \J\x\<ri J \J\x\<r) 

and using that (j) q ^ n with n > are uniformly bounded (by (tt — l) -1 / 2 ) and that 



dx 4> q ,ni x ) 



dx \ip(x,q,t)\ 



/ dq R(r], q) < Cr] 
Jr 



for a suitable constant C. Similarly, using Cauchy-Schwarz and J dx<pg n (x) = 1, 

Us M <l*,l „«)!', 

Jr Jr Jr)<\x\<n 

so, thanks to Proposition [31 



limsup / dq\Q n (q,t)\ 2 < Cr] , 

t->oc JR 



and the claim follows because rj > is arbitrary. □ 

The above argument fails if n = 0, because 9) o is n °t uniformly bounded in q < 0. 



(83) 



(84) 



(85) 



1. 



G. About the Band potential. 



Here the 3d term in the band potential (J39J) is estimated. A standard perturbative calculation 
yields: 



1nl(q) = I dx 4>q,n( X ) 



d<Pq,l(x) 

dq 



a W l {q)-W n {q)' {nti) 



(86) 
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and j nn = 0; so, thanks to orthonormality and completeness of {4> q , n }: 

H^)' = (ST, 

It is easy to see that \4>q, n {0)\ < (vr - 1)~ 1/2 whenever n > 0, that W n {q) > \r? , and that 
1 09,0 1 ~ V~ a 1 f° r 9 — ^ ~~ 00 • Using this and the asymptotic form of Wo(q) given in ([35]) . ((871) is 
found to be 0(q~ 2 ) for g — > — oo and ~const. as q — > +oo. 



H. Proof of Proposition [5} 

The following notations will be used. For sufficiently small rj > 0, Aj„ = Jj \ Dj where 
Dj C Jj is an arc of size r] > centered at Oj\ and = UjAj^. P v will denote the projector 
of L 2 (§) ® L 2 ^) onto L 2 (A T? ) ® L 2 (M Ar ), and P ii7? will denote the projector of L 2 (J i ) L 2 ^) 
onto L 2 (Aj :V ) L 2 (R N ). 
Existence of fi ± = n±(H^>,'H) entails that: 

lim ( e - m V .-Pfi e~ iW V) = I™ (e- wC&) *0 ± ^ , P„ e^^Oi^) 

i— >±oo f^±oo 
TV 

= MoX {e ~ m{b)t ^ ' P ^ e-** W *n^) ■ (89) 
^ °°i=i 

The quantity of which the i — >• oo limit is taken in the above equations is the probability of finding 
the particle in A v at time t. Each subspace Pj($j) is invariant under the evolution ruled by H^', 
so the sum in the last line is equal to 

2V 

Y,{e- mib)t P^ ,Pj,r, e- mib)t P^) . (90) 

3=1 

Each term in the sum is a probability of finding the particle at time t in Jj^ with the evolution 
e —iH b t j n invariant subspace this evolution is that of a single-oscillator model, so, thanks 
to Proposition [3j each term in the sum tends to as t — > oo. Hence, so does the probability of 
finding the particle in A„ for all r\ > 0, which is equivalent to the thesis. □ 
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